Quantum Shuttle in Phase Space 
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We present a quantum theory of the shuttle instability in electronic transport through a nanos- 
tructure with a mechanical degree of freedom. A phase space formulation in terms of the Wigner 
function allows us to identify a cross-over from the tunnelling to the shuttling regime, thus extending 
the previously found classical results to the quantum domain. Further, a new dynamical regime is 
discovered, where the shuttling is driven exclusively by the quantum noise. 
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Advances in microfabrication technology are push- 
ing today's microelectromechanical systems (MEMS) to- 
wards the nanometer regime, and the emerging new tech- 
nology of nanoelectromechanical systems (NEMS) is ex- 
pected to play an important role in the future. The ubiq- 
uitous quantum mechanical effects affecting the perfor- 
mance of these devices present many theoretical chal- 
lenges, only few of which have so far been addressed in 
the literature. The purpose of this Letter is to present a 
fully quantum theory for an electromechanical instabil- 
ity in a generic NEMS device,the single-electron shuttle 
first studied by Gorelik et al . 

This device consists of a movable single-electron tran- 
sistor (SET) and exhibits an electromechanical instability 
from the standard tunnelling regime to a new regime in 
which the SET oscillates and carries an integer number 
of electrons per a cycle (shuttle regime). Since the orig- 
inal suggestion 0, there has been increasin g in terest in 
the shuttle phenomenon 0, @, 0, 0, II S H1I3 : e.g. by 
incorporating the shot noise due to the electron transfer 
, gate effects , the coherence effects in the electronic 
subsystem and strong dissipation of the oscillator en- 
ergy 0. Only very recently the quantum mechanical 
treatment of the oscillations in various modifications of 
the shuttle setup has been considered . 

The classical theory of shuttle transport has been used 
00 to describe the experiments on C 60 single-electron 
transistor , where the oscillations of the center of mass 
of the molecule were found to be important. However, 
also explanations based on incoherent phonon assisted 
tunnelling theory which do not take into account the cor- 
relation between the coherent oscillator motion and the 
electron transfer seem to yield reasonable predictions for 
the I-V curves [HEl. 

Therefore, a fully coherent quantum mechanical treat- 
ment of the oscillator which does take into account pos- 
sible correlations and which would, in principle, allow 
to join the two approaches into a unified framework 
is a most desirable theoretical task. In this study we 
present an attempt for such a treatment for a simplest 



model supporting the shuttling transition in the classical 
regime. A complementary study to ours already exists 
I Id) which uses a similar model, albeit without any me- 
chanical damping. However, the quasiclassical expansion 
of the tunneling term used in that work does not give ac- 
cess to the purely quantum phenomena discussed below. 

We demonstrate that the characteristic strong correla- 
tion between the oscillator motion and the electron trans- 
fer persists even in the quantum regime. The noise gen- 
erated by various sources (shot and thermal, both having 
quantum components) is found to be very important for 
the phenomenon. Not only does it smear the classical 
transition found in jl| into a crossover with a consid- 
erably shifted position in parameter space compared to 
the classical "mean field" study 0, but the quantum 
component of the shot noise also generates the shuttle 
instability even in the classically stable region with zero 
electric field. 

Using the generalized master equation approach sug- 
gested in we study a simple model motivated by sev- 
eral previous studies 0, |j, [ljj, [ijj • Namely, we consider 
an oscillating nanoscopic grain with only one electronic 
level (strict Coulomb blockade regime) coupled to two 
leads. The oscillator degree of freedom is treated fully 
quantum mechanically in our approach. We also account 
for the oscillator damping, the thermal noise due to the 
oscillator-bath coupling, and the shot noise due to the 
quantized electron transfer. 

The Hamiltonian of the model reads 

H = H osc + EqcIcq + ^ ( e ka ~ Ha)c\ a C ka + H B 
k;a=L,R 

-eExc' c + ^ {Tkao{x)c\ a C(j + h.c.) + H mt 

k;a=L,R 

(1) 

where the first three terms describe the free evolution 



of a linear harmonic oscillator H n 



+ 



the 
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single electronic level on the grain, and the two nonin- 
teracting leads symmetrically biased by voltage —V (i.e 
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[Xl = ^j-, Pr = — 3" with e, V > 0; all electronic ener- 
gies are measured with respect to the equilibrium chem- 
ical potential of the leads), respectively. The term Hb 
describes a generic Ohmic heat bath |l4| . The first in- 
teraction term accounts for the electrostatic coupling be- 
tween the position of the grain and the occupation of the 
grain electronic level; E is the electric field caused by 
the bias applied between the leads and/or gate voltage 
0, E3 • This electric field causes a shift of the charged 
oscillator equilibrium by d — The second term is 

the coupling of the leads to the single electronic state on 
the grain with oscillator-dependent hopping amplitudes 
Tkw(x) = t L exp(-f), T kR0 (x) = t R exp(f) where A is 
the electron tunnelling length. The last term H[ nt de- 
scribes the coupling of the heat bath to the oscillator 
which is linear in the oscillator coordinate |l4|. The net 
effect of this coupling is the oscillator mechanical damp- 
ing characterized by a constant 7 and the Langevin ran- 
dom force that the bath exerts on the oscillator and which 
depends on the temperature of the bath T. 

By projecting out the leads and the thermal bath 
we arrive at a Markovian generalized master equation 
(GME) for the density matrix p(t) of the system com- 
posed of the grain and the harmonic oscillator. Since the 
electronic transfer rate in the shuttling regime is compa- 
rable with the frequency of the oscillations we employed 
the singular coupling limit 0,0] appropriate for rapidly 
decaying systems. This approximation is valid for bias 
much larger than any energy scale of the system, i.e. 
eV » nui,e . Moreover, we also assume eV » fc^T 
which is reasonable in the present physical context. The 
mechanical damping due to the heat bath is on the other 
hand treated within the standard weak coupling theory 

Our GME reads 



p{t) = Cp{t) = (£ coh + £driv + £damp)p(i) 



(2) 



with 



CcohP = ~f\ H osc + eocJc - eExclc , p] , (3) 
£drivP = --^-(cocle'^p - 2cje~* pe'^c + pe~^ Cf)C t) 



_ -TT ( c o c o e 2 " P ~ 2c oe * pe * cl + pe ^ cJcq) , 



(N + l/2)[x,[x,p}} 



(4) 



(5) 



Here, [,] and {,} denote the commutator and the anti- 
commutator, respectively and N = (exp(-^y) — is 
the mean oscillator occupation at the bath temperature. 

The first term in Eq. (J2J) is the free coherent evo- 
lution of the oscillator and the grain level while the 
second one describes the transfer of electrons via the 



oscillator-position-dependent tunnel junctions from the 
left reservoir to the grain level (the term with T^) and 
from the level to the right reservoir (the term with Tr). 
The opposite processes can be neglected due to the as- 
sumption eV 3> huj,€Q,kBT. The transfer rates equal 
Fl.r — ^\tL,rt\ 2 T^L,R with the constant densities of 
states of the leads T>l,r- The third term accounts for the 
interaction of the oscillator with the heat bath. Trans- 
lational invariance of the damping, positivity of the den- 
sity matrix, and relaxation towards canonical equilibrium 
cannot be achieved simultaneously with any Markovian 
damping kernel and one has to sacrifice at least one of 
these properties. The most physical choice is to relax the 
positivity 0] . We checked the magnitude of breaking of 
the positivity in our calculations and found it irrelevant. 
Moreover, it only occurs for large values of 7 out of the 
shuttling regime. 

It can be shown that the electronic off-diagonal el- 
ements of the density matrix are decoupled from the 
diagonal ones and, moreover, decay to zero in the sta- 
tionary state. Therefore, it is sufficient to consider only 
the electronic diagonal elements: Poo(*) = (0 |p(^)|0) and 
p n (t) = (l\p(t)\l), where |1) = cj|0>. These objects 
are still full density matrices in the oscillator space and 
satisfy 

Poo(t) = ^r[H OBC ,p o(t)} - ^(e~^ p 00 {t) + poo(t)e~^) 
in 2 

+ T R e%pti(t)e% + £damp Poo(t) , 

1 _* 

Pn(t) = —[H osc - eEx,pn(t)]+T L e x p o{t)e * 
in 

— (e"*"pn(i) + pn(t)e~) + £dam P pu(t). 

(6) 

From the continuity equation for the electronic charge 
we may deduce the following formula for the stationary 
current through the grain (flowing from the left to the 
right lead): 

/s tat = eTiT^fe-^pg,'*) = eT R Ti osc (e^ pff) . (7) 
The trace is carried out over the oscillator basis and 

Pnn = Km^ooPnn(t) ■ 

We solved numerically the stationary version of the 

above equations 10: = Yll=i Lkipf at where the 
column vector pf at consists of the matrix elements of 
/5gg at , p\ * at and the (super-)matrix Lki of the dimension 
2N 2 x 27V 2 contains the appropriate coefficients of the 
linear system (JBj. The density matrix was represented 
in the harmonic oscillator basis which was truncated by 
taking up to N = 100 lowest states which yields satisfac- 
tory numerical convergence in our parameters range. We 
determined the unique (for T R = Tl = T ^ 0) null vec- 
tor of the supermatrix L by using the Arnoldi iteration 

m 
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T= 0.25 y= 0.1 y= 0.05 




FIG. 1: Phase space picture of the tunnelling-to-shuttling crossover. The respective rows show the Wigner distribution functions 
for the discharged (Woo), charged (Wn), and both (Wtot) states of the oscillator in the phase space (horizontal axis - coordinate 
in units of xq — \J h/muj, vertical axis - momentum in h/xo). The values of the parameters are: A — xq, T — 0, d = 0.5xo, Y — 
0.05hui. The values of 7 are in units of hui. The Wigner functions are normalized within each column. 



Our approach differs from the one used in |fj where the 
explicit time integration scheme of the time-dependent 
equation analogous to © was used to determine the sta- 
tionary density matrix. We applied our method to the 
model of and recovered the I-V curve results reported 
there. Using the phase space analysis detailed below we 
fully confirmed the shuttling interpretation based on the 
indirect evidence from changes of the I-V curves with 
changing parameters. 

The I-V curve (or other dependencies of the current on 
some parameter) alone yields only an indirect evidence 
of shuttling and may actually not be decisive whether 
the system is shuttling or not (see e.g. 01 versus 
Therefore, it is preferable to consider quantities which 
depend also on the state of the oscillator. An excellent 
visualization tool 0] for the description of the joint elec- 
tronic and oscillator properties are the Wigner functions 
(n = 0, 1) 

(8) 

yielding the charge-resolved quasiprobability distribu- 
tions of the oscillator in the phase space. These functions 
provide us with a clear evidence of the transition from 
the incoherent tunnelling to the coherent quasiclassical 
shuttling behavior with decreasing damping coefficient. 



Our focus is in the quantum effects on the shuttling 
transition and, in particular, whether there is any transi- 
tion in the quantum regime at all. Therefore, we work in 
the strictly quantum regime where the tunnelling length 
is comparable to the zero uncertainty of the oscillator, 
A ~ xq = yjh/muj, and where only a relatively small 
number of oscillator states is excited. 

In Fig.Q]we depict the Wigner functions Woo, Wn, and 
Wtot — Woo + Wu showing the crossover from the tun- 
nelling to the shuttling regime with decreasing damping. 
In the tunnelling regime (large 7) the oscillator is located 
around the origin (or shifted origin when charged) with 
no particular correlation between its charge state and 
momentum (Wigner functions are centered around the 
origin with some "quantum fuzziness"). This is consis- 
tent with the quantum incoherent tunnelling picture. On 
the other hand, in the shuttling regime (small 7) the os- 
cillator orbits almost classically (ring-like shape of W to t 
with a hole around the origin), and shuttles the charge 
on its way from the left to the right lead and returns 
empty back (half-moon shapes of Woo, Wu). The corre- 
lation between the charge state and the mechanical mo- 
tion is very strong. In the crossover region (medium 7) 
we can see that both regimes of transport are contribut- 
ing additively (ring-like shape plus an incoherent peak 
around the origin of Wtot)- The classical "mean field" 
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FIG. 2: 1 — 7 curve. The 7-dependence of the stationary 
current through the grain for different transfer rates and elec- 
tric fields. Their values are d = 0.5xo,r = 0.05/kj (pluses; 
corresponds to Fig. 0, d = 0.5a;o,r = O.Olhuj (circles), 
d = 0.0, T = 0.05ftw (asterisks), d = 0.0, T = 0.01ft (crosses). 
Other parameters are A = xo,T = 0. The current is in units 
of euj while 7 in Hu. 

sharp transition of between the tunnelling and the 
shuttling regime is smeared into a crossover due to the 
noise. Also the position of the crossover is substantially 
shifted with respect to the classical values (7 cr oss is ~ 5 
times larger than the classical value) which is attributed 
to the deeply quantum (and therefore noisy) regime. The 
classical picture is expected to emerge in the quasiclassi- 
cal limit rf,A>io. 

In Fig. |21 we plot the 7-dependence of the stationary 
current through the grain for different transfer rates T 
and electric fields (measured in d). For the two curves 
with d 7^ we can see the rise of the current in the 
crossover region from the tunnelling-limited values pro- 
portional to T to shuttling-mediated quantized value of 
one shuttled electron per each cycle (1/27T ~ 0.16 in our 
units; independent of T) in agreement with the classical 
results. More surprisingly, also the results for the case 
of zero electric field d = show clear signs of shuttling 
crossover in the I — 7 curve for small enough mechanical 
damping. Classically there is no shuttling transition for 
d = regardless of the other parameters values and 
the 7—7 curve is constant. Therefore, this shuttling must 
be driven purely by the quantum component of the shot 
noise (proportional to h). Also the phase space pictures 
of this regime reveal onset of shuttling transport. 

Finally, we comment on the effect of the temperature. 
For a nonzero temperature the shuttling transition within 
our model is facilitated by the increase of the mechanical 
noise driving the transition. The deteriorating effect of 
the temperature on the transition is not included in the 
model due to the high bias assumption. The development 
of the theory for a finite bias is under way. 

To summarize, we have presented a quantum theory 



of the shuttling transition in fully developed Coulomb 
blockade regime. Using the Wigner functions as a phase 
space visualization method we have exhibited a clear 
crossover from the tunnelling to the shuttling regime of 
the transport as a function of the mechanical damping 
parameter. The effect of noise on the transition in the 
deeply quantum regime (A ~ xq — ^Jh/muj) is pro- 
nounced and can even trigger the transition in the clas- 
sically stable regime with the zero electric field. 
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